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The semiclassical quantization of cyclotron orbits for two-dimensional Bloch electrons in a coupled 
two band model with a particle-hole symmetric spectrum is considered. As concrete examples, we 
study graphene (both mono and bilayer) and boron nitride. The main focus is on wave effects - 
such as Berry phase and Maslov index - occurring at order h in the semiclassical quantization and 
producing non-trivial shifts in the resulting Landau levels. Specifically, we show that the index shift 
appearing in the Landau levels is related to a topological part of the Berry phase - which is basically 
a winding number of the direction of the pseudo-spin 1/2 associated to the coupled bands - acquired 
by an electron during a cyclotron orbit and not to the complete Berry phase, as commonly stated. 
As a consequence, the Landau levels of a coupled band insulator are shifted as compared to a usual 
band insulator. We also study in detail the Berry curvature in the whole Brillouin zone on a specific 
example (boron nitride) and show that its computation requires care in defining the "k-dependent 
Hamiltonian" H(k), where k is the Bloch wavevector. 

I. INTRODUCTION 

The dispersion relation of Bloch electrons in two dimensional (2D) crystals generally exhibit regions of closed 
orbits in reciprocal space. As a consequence, it is expected that applying a perpendicular magnetic field gives rise 
to quantized cyclotron orbits and the corresponding Landau levels. A semiclassical approach to obtain these Landau 
levels consists of first computing the area of the classical cyclotron orbits and then imposing the Bohr-Sommerfeld 
quantization condition in the form suggested by Onsager for Bloch electronsi. The semiclassical quantization condition 
(see Appendix E| for a cyclotron orbit C reads: 

5(C)/| = 2^[n + 7] (1) 

vifhere 5(C) = // (Pk is the fc-space area enclosed by the cyclotron orbit, k is the (gauge-invariant) Bloch wavevector, 
Ib = y^n/eB is the magnetic length, — e is the electron charge and n is an integer. The quantity 7 is called a phase 
mismatch (0 < 7 < 1) and is not given by the semiclassical quantization rulei. The precise determination of 7 
requires the inclusion of wave effects and therefore to include terms of order h in the semiclassical expansion. For free 
electrons, and for a single (uncoupled) band of Bloch electrons, 7 = 1/2 as a result of the presence of two caustics on 
the cyclotron orbili^i^. The number of caustics on an orbit is known as the Maslov index— 1^. From the dependence of 
the cyclotron surface S{C) on the energy e, one can usually rewrite the above quantization condition as: 

5(£)/| = 27r[n + 7i] . (2) 

Then by inverting S'(e), one obtains the (semiclassical) Landau levels 

2tt 

e„ = S-^[-^{n + jl)] = function[S(n + 7^)] (3) 

where n is now interpreted as the Landau index. Usually, the shift 7^ is trivially equal to the phase mismatch 7 
introduced above. For example, the Landau levels for a free electron of mass m and dispersion relation e — fi^k^ /2m 
are given by a harmonic oscillator e„ = (n -I- l/2)heB/m, and 7^ = 7 = 1/2 in that case. Indeed the area of the 
cyclotron orbit is S{C) ~ Trfc^ and therefore iS'(e) = 2'K'melY? such that S[e)l\ = 27r(n -I- 1/2). One of the goal of this 
paper, it to show that these two quantities, 7 and 7^, are not necessary equal. 

A relation between the phase mismatch 7 and the nature of the electronic Bloch functions was obtained by Rotb^. 
She found that 7 can depend on the cyclotron orbit C and that 7(C) can be related to a quantity r(C) later identified 
by Wilkinson^ as a Berry phased acquired by the Bloch electron during a cyclotron orbit C, see also Ref. 0. The 
relation reads 

7(C)=7M + 7B = ^-^^ (4) 

Z ZTT 
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where 7jv/ = 1/2 refers to the Maslov index contribution and 7b = —T{C)/2n to the Berry phase contribution. The 
Berry phase is given by 

T{C)^i j dk- {uk\VkUk) (5) 

in terms of the Bloch function Uk{r), where k is the gauge- invariant Bloch wavevector, and is computed along the 
cyclotron orbit C. 

A case which is particularly interesting from the perspective of semiclassical quantization is that of coupled bands. 
In the present paper, we will restrict to two coupled bands with electron-hole symmetry, having in mind the examples 
of graphene (the two bands touch at two inequivalent valleys known as Dirac points) and boron nitride (the two 
bands are separated by a gap). In a two coupled band system, the Bloch electron is endowed with a pseudo-spin 1/2 
associated with the freedom of being in the two bands and its wavefunction is therefore a bi-spinor. In the context of 
graphene or boron nitride, this internal degree of freedom is usual called "sublattice pseudo-spin" as it results from 
having two inequivalent sites A and B in the unit cell. 

The Landau levels of electrons in graphene were first obtained by McCluro^, who performed a fully quantum 
mechanical calculation and obtained the now well-known behavior 

£n,a = C(vV 2neBh (6) 

where a = ±1 is the band index and v is the constant Fermi velocity. From a semiclassical perspective^, this result 
and the value 7l = that it implies - via equation ([3|) - seem to imply that the phase mismatch is now 7 = instead 
of the usual 7 = 1/2. Using the Roth- Wilkinson relation (|4]), Mikitik and Sharlaii were able to show that this value 
of 7 can be attributed to a Berry phase F = tt, which exactly cancels the Maslov contribution. This Berry phase of 
TT is actually the phase appearing in Hilbert space when rotating a bi-spinor by an angle of 27r in real space, which is 
familiar in the context of spin 1/2 physics, see e.g. Ref. [^. The conclusion is that 7^ = 7 = in graphene and this 
value has indeed been observed in Shubnikov-de Haas and quantum Hall effect measurements^^. 

In this paper, we focus on a situation where the assumption that 7l = 7 gives a wrong result for the Landau levels. 
This is, for example, the case of boron nitride^, whose low energy effective theory is that of massive Dirac fermionsH, 
with a dispersion relation e{k) = a-^/A^ -I- {vhkY where the gap 2A = Eb — £n is the energy difference between a 
boron and a nitrogen 2pz atomic orbital. Haldan&i^ computed the Landau levels of two dimensional massive Dirac 
fermions quantum mechanically and found^: 

£n.a = a + 2neBhv^ . (7) 

From a semiclassical perspective Haldane's result raises the following question. When comparing masslcss 
(graphene) and massive (boron nitride) Dirac fermions, it appears that in both cases, the Landau level shift is 
the same 7^ = . However, as we will show below, the Berry phase F(C) depends on the magnitude of the gap and 
is therefore different in both cases. We are led to conclude that the two quantities 7 and 71, are different in this case. 
The main goal of this paper is to relate these two quantities. We will show that while 7 entering the quantization 
dl]) of cyclotron orbits is correctly related to the Berry phase, jl entering the energy quantization ^ is related to 
a topological part of the Berry phase, which is essentially a winding number of the pseudo-spin 1/2. The key point 
to understand this difference is to account for the orbital magnetization of Bloch electrons and we will see that when 
quantizing the cyclotron orbit one needs to consider the change in energy due to this magnetization. In the Onsager 
quantization condition, the contribution of the orbital magnetization exactly cancels the non-topological part of the 
Berry phase. As a consequence, the Landau index shift 7^ is only the topological part of the Berry phase and not 
the whole Berry phase. 

Recently Carmicr and UUmoi^ have computed the semiclassical Green functions for similar systems (graphene, boron 
nitride, etc.). With a different approach, they reach essentially the same conclusion as ours: the phase appearing in 
the Landau levels is in gene ral not the complete Berry phase but is what they call the semiclassical phase. In the 
same vein, see also Ref. lla. Another related work, which appeared recently, is that of Gosselin, Bcrard, Mohrbach 
and Ghosbii. Compared to these two works, on the one hand, ours treats more general coupled band systems and do 
not require to take the continuum limit but also applies to discrete models on a lattice. For example, the quantization 
of cyclotron orbits can be performed anywhere in the Brillouin zone and not only close to specific points. On the 
other hand, and contrary to Ref. [iB-ES we restrict ourselves to homogeneous systems. 

The structure of the paper is as follows. In section II, we review the semiclassical description of Bloch electrons at 
order H having in mind the quantization of cyclotron orbits and introduce various Berry quantities. Section HI is the 
core of the article. It contains a study of a two coupled band model for which we show that the Landau index shift 
7i is related to a winding number and not to the complete Berry phase. Then in the following sections, we consider 
several examples: a tight-binding model for boron nitride (IV), massive Dirac electrons (V), massless Dirac electrons 
(VI), and eventually chiral electrons of bilayer graphene (VII). The conclusion is presented in section VIII. 
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II. SEMICLASSICAL DESCRIPTION OF A BLOCK ELECTRON ON A CYCLOTRON ORBIT 

111 this section, we review known results about the semiclassical description (including terms at order h) of Bloch 
electrons in a crystal under the influence of a magnetic field. When describing a Bloch electron confined to a single 
band, the presence of other bands shows up in the semiclassical equation of motions at order h (classical order being 
hP) in the form of Berry phase type corrections. Our goal is to discuss the effect of these corrections on the quantization 
of cyclotron orbits. It should be kept in mind in the following that the electron is described by a wavefunction which 
is a bi-spinor - because of the band structure - and that the true spin is neglected. For a general review see Ref. fisl . 

A. Semiclassical equations of motion for a Bloch electron in a magnetic field 

One wayi^ii^ of obtaining the semiclassical equations of motion for a Bloch electron in a uniform magnetic field B 
is to study the motion of a (typically Gaussian) wavepacket of Bloch waves restricted to a single band (indexed by a) 
of average position rdt), average crystal momentum hqj(t) - is the average Bloch wavevector - and fixed width. 
The width of the wavepacket should be larger than the lattice spacing and much smaller than the typical length scale 
on which the external fields (e.g. magnetic and electric) vary. One then uses the time-dependent variational principle 
to obtain an effective Lagrangian for the independent variables Tc and q^. Minimizing the action with respect to these 
variational parameters one obtains the following equations of motion: 

hkc — —ere X B (8) 

and 

rc = hr^Vk^ea-kc-x^a{kc) (9) 

where hkc ^ hq^ + eA(rc) is the average gauge-invariant crystal momentum^, A is the vector potential and — e < 
is the electron charge. The Berry curvature ^a{kc) is defined below. The electron energy is 

ea(fcc)=ea,o(fcc)-A^a(fcc)-B (10) 

where ea.o{kc) is the band energy in absence of a magnetic field and AA.a{kc) is the orbital magnetic moment of the 
Bloch electron (also defined below). 

Compared to the usual equations of motion of Bloch and Peierls [see e.g. Ref. [l^l obtained at order W , there are 
two additional terms in Eq. (|8|9p . which appear at order h. One is the so-called anomalous velocity —kc x ^a{kc) ■ 
It is a kind of Lorcntz magnetic force but in fc-space and due to Berry curvature JIq, rather than to a real magnetic 
field. It takes into account the effect on the average velocity of virtual transitions to other bands a' ^ a. The other 
is the magnetization correction to the band energy, which gives the energy of a Bloch electron in a magnetic field 
as Ea = £a,o ~ AIq, ■ B . The correction to the band energy is the extra magnetic energy due to the coupling of 
the orbital magnetic moment Ala(fcc) to the external magnetic field. This orbital magnetic moment comes from the 
self-rotation^ of the wavepacket and exists because of the finite wavepacket width, which cannot be made arbitrary 
small due to the restriction to a single bandit. Just as the Berry curvature, the orbital magnetic moment is also an 
effect of virtual transitions to other bands. In the present work, we neglect the electron spin and therefore do not 
discuss the Zeeman effect. Note that the magnetization appears as resulting from an internal structure, which is not 
the electron spin but rather the pseudo-spin related to the two coupled bands. 

Berry-type corrections, such as fjQ,(fec) and A^a(fcc), appear at first order in the external field. Therefore, to this 
order, ^a{kc) ~ Q.a{qc) and J^a{kc) ~ J^aiQc) and it is therefore not important to distinguish between k and q 
when computing these quantities. 

B. Berry-ology 

Here we consider the effective dynamics of an electron restricted to a single band and define several quantities 
related to a Berry phase appearing because of the coupling between bands. In particular, we consider a 2D crystalline 
system described by a Hamiltonian H containing only two bands (band index a = ±1). Typically, we think of a 
tight-binding model with two sites in the unit cell. In the Bloch basis it reads: 



H = ^ea,o{k)\k,a){k,a\ 



(11) 
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where ea.o{k) is the band energy, \k,a) = exp(ifc • r)\uk,a) is a Bloch state and r is the complete position operator 
(and not just the Bravais lattice position, e.g.). Its wavefunction is (pk.aif) = (^l^jQ^) = exp(ifc • r)uk^a{f): where 
Uk,a{'r') is the Bloch function. In the case of two bands, <fk,a{''') and Ufc,a(^) are bi-spinors (in sublattice space). Next, 
we perform a unitary transform to define a fc-dependent HamiltoniaiJ^: 

H{k) = cxp{—ik ■ r)H cxp{ik ■ r) (12) 

The wavcvector fc is a parameter spanning the first Brillouin zone and on which the Hamiltonian Hik) depends. By 
virtue of the unitary transform exp(— ifc • f), one has H{k)\uk,a) — £afl{k)\uk^a) ■ Using the projection operators 
P{^) = \uk.a){uk.a\i we also define the following 2x2 fe-dependent Hamiltonian: 

H{k) = P{k)H{k)P{k) - ^eo^a(fc)|Mfc,a)(«fc,o| (13) 

a 

which is the restriction of H{k) to the k subspace. For more details on the three different types of Hamiltonians we 
are using [i?, H{k) and H{k)] see Appendix IB] 

Following the general result of Ref.— , the Berry phase acquired by a Bloch electron on a cyclotron orbit^i^ C in 
the band a is: 

Ta{C) = <j> dk ■ i{uk,a\^kUk,a) (14) 

Note that in general this quantity depends on the cyclotron orbit C . The Berry connection (equivalent to a fc-space 
vector potential) in the band a is given by: 

Aa{k) = i{uk,a\V kUk,a) (15) 

SO that the Berry phase appears as an Aharonov-Bohm phase in fe-space. The corresponding Berry curvature (equiv- 
alent to a fc-space magnetic field) is 

f^a(fc) = Vfe X ^„ = n^e^ (16) 

where 

^a{k) = dk^Ay - dkyAx = i[{dk^u\dkyU) - {dkyU\dk,u)] (17) 

It can also be written as: 

Cta{k) = i{VkUk,a\ X \VkUk,a) (18) 

Another useful formulation, especially convenient when performing numerical calculations as, contrary to Eq. ()18p . 
it docs not require the Bloch wavefunctions to be single- valued in parameter spaced, is: 

O /7^ ■ {uk.a\dk^H{k)\uk^a'){uk.a'\dkyH{k)\uk,a) , 

f^afe =«V r TT-. ; TTTil + C.C. 19 

^ £a,0(fe) - £a',0(fe)r 

It shows explicitly, that the Berry curvature is due to the restriction to a single band a and to the resulting virtual 
transitions to other bands a' ^ a. 

The orbital magnetic moment Ala = M-a^z of a Bloch electron described by a wavepacket of average position rc 
and average gauge-invariant crystal momentum hkc restricted to the band a isi^ii^: 

6 C 

Maikc) = -T^{{r - Vc) X p) = -i — (Vfe^Ufc^.a! X [e„^o - H{kc)]\V k,Uk^,a) (20) 
2m 2h 

where the average in the first expression is taken over the wavepacket, p is the canonical momentum operator, and m 
is the bare electron mass. As the Berry curvature, this quantity also has an expression revealing the virtual transitions 
to other bands: 

■ e {uk,a\dk^H{k)\uk,a'){uk,a'\dkyH{k)\uk.a) , 
•^a{k)^t— > TTT TTV '■ ^C.C. (21) 
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This shows that in the case of a single isolated band, both the Berry curvature and the orbital magnetic moment 
vanish. Note that both quantities depend on the off-diagonal (in band index) matrix elements {uk a\^ kH {k)\uk,a') 
of the velocity operator. 

In the particular case of a two-band model with electron-hole symmetry, the orbital magnetic moment is directly 
related to the Berry curvature: 

Mo, = |ea,of^a (22) 

This relation was already obtained in Ref. [12 and we present a proof in Appendix [C] 

According to general symmetry arguments^^, the Berry phase and the magnetization of a single band should vanish 
in a crystal which is inversion and time reversal invariant. Indeed, time-reversal symmetry implies 11(— fc) = — ri(fc) 
and inversion symmetry implies ft{—k) = r2(fc). 

All the above definitions are valid for an electron in a single Bloch band, which is well separated from other bands. 
We will nevertheless apply them in the case of touching bands (such as graphene at its Dirac points) remembering 
that the correct procedure is to calculate these quantities in presence of a finite gap A and to send it to zero at the 
end. 



C. Cyclotron orbit, phase mismatch and Landau index shift 

In the following, the aim is to quantize the cyclotron motion in order to find the Landau levels. Classically, a free 
electron in a uniform and constant magnetic field performs a motion at constant energy in a plane perpendicular to 
the magnetic field. For a Bloch electron, the classical cyclotron orbit is a cut at constant energy in the band structure, 
i.e. an iso-energy line ea{k) = constant. The semiclassical quantization of a cyclotron orbit is explained in detail in 
the introduction - see equations ([1]), and ([5]) - we therefore do not recall it here. Nevertheless, we would like to 
precise the definition of the Landau index shift 7l, which is related, but not identical, to the phase mismatch 7(C) 
appearing in the Onsager semiclassical quantization condition ([T]). The Landau index shift appears in the energy 
quantization condition ([2]). It can also be defined via the exact Landau levels e„ by taking the semiclassical limit 
(n ^ 1, keeping terms of order n and n*^): 

En ~ function[_B(7i + 7l)] (23) 

where n corresponds to the dominant term, of order and to the first correction, of order n'^ ^ To be 

more precise, imagine expanding the exact Landau levels as a decreasing series in powers of n: e„ = ao?i' + ai?i'^^ + 
a2'n}~'^ -|- . . .. Keeping only the two first terms in the semiclassical limit n 3> 1, one obtains £„ w avi[n^ +ain'^^/ao] ~ 
ao[ri + (111 (floO]' = (io[n + 7l]', which defines the Landau index shift 7^ = ai/(ao/) modulo 1. 

Often both quantities 7(C) and 7^ are equal and are usually not distinguished. The insight here comes from 
recognizing that both quantities can be different as 7(C) may depend on the precise cyclotron orbit, whereas 71, is a 
constant. 



III. SEMICLASSICAL QUANTIZATION OF CYCLOTRON ORBITS IN A COUPLED TWO-BAND 

MODEL 



In the following, we perform the semiclassical quantization of the cyclotron orbit for a Bloch electron in a two-band 
model and obtain the relation between 7(C) and 71,. We consider a coupled two-band Hamiltonian with a particle- hole 
symmetric spectrum H . As explained in the previous section, we then perform a unitary transform exp(— jfc • f) to 
obtain a parameter-dependent Hamiltonian H(k) and then project on the k subspace to obtain a 2 x 2 Hamiltonian 
[in the following ft = 1] : 

^(^) - ( r%) -A ) (24) 

where k is the Bloch wavevector in the first Brillouin zone (BZ). The function f{k) is usually obtained as a sum over 
hopping amplitudes in a tight binding description. Time-reversal symmetry imposes H{—k)* = H{k) and therefore 
f{—k)* = f{k). Note that Bloch's theorem imposes that |/(fc + G)\ = |/(fc)| for any reciprocal lattice vector G. 
However it does not require that f{k + G) = f{k). An important assumption here is that the diagonal term A does 
not depend on the wavevector and can therefore be interpreted simply as an on-site energy. This term explicitly 
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breaks the inversion symmetry. Introducing the energy spectrum eo{k) = a^/ + |/(fc)p, where a = ±1 is the band 
index, and the azimuthal P{k) and polar 9{k) angles on the Bloch sphere, such that cos/3 = A/|eo|, sin/3 = |/|/|eo| 
and 9 = — Arg/, the Hamiltonian can be rewritten as 

j^,... . .f cos/3 sin/3e"*^ \ 

H{k) = \eo\ ■a^^e « (25) 



sin /3e — cos /3 

The eigenfunction of energy eq = a|eo| is i'ir) = Ufc(r)e*'° '' where the Bloch spinor is 

, , f cos(/3/2) \ 

1"'=-^ = i sin(/3/2)e''' 1 if a = +1 



-sin(/3/2)e-''' 
cos(/3/2) 



ifa==-l (26) 



The Berry connection is given by 



and the corresponding curvature is 



^= -asin^ |Vfc6' (27) 



^='^VkCosPxVke = -'^smp{dkJdk^e-duJdkJ)e,. (28) 

An important simplification occurs in the calculation of the Berry phase F because the cyclotron orbit C is travelled 
at constant energy and the diagonal term A is independent of the wavevector. As a consequence, the azimuthal angle 
/3 is a constant along the trajectory. Indeed cos/3 = A/|£o| and sin/3 = y's^ — A^/jeol are both functions of £o only. 
Therefore the calculation of the Berry phase along a cyclotron orbit is easily performed: 

r(C) ^ <j> dk • A = -asin^ ^ j '^^ ' ^'^^ " t^Wc[1 - cos/3] (29) 

where Wc = §^d9 /2t: is the winding number, which is a topological invariant. Indeed the relevant mapping is 
from a cyclotron orbit in the Brillouin zone to a circle (because /3 is fixed) on the Bloch sphere: therefore, the relevant 
homotopy group is 7ri(S'^) = Z. Note that (i(|eo|r)/rf|£o| = constant = ttWc- We call this quantity the topological 
Berry phase. It is a local quantity as it depends on the precise path C. The winding number Wc counts the total 
charge of the vortices in 0, which are encircled by the cyclotron orbit (see Figure 2). Note that this topological Berry 
phase is not directly related to the Chern number, which is the Berry curvature integrated over the entire BZ22. 
Starting from the Onsager-Roth relation (see Eq. (IlKISp *) 

5(£o)^| = 27r[7i + i] - r(C) (30) 

where Eq is the band energy in zero magnetic field, we search the quantization of S{e) where e is the energy in presence 
of a magnetic field. Using the relation between the energy and the curvature Eq = e + M.B with = ceq^, we obtain 

S{e,)ll^S{e)ll + n[eo)\e,\^^ (31) 

In the previous equation, we introduced the Berry curvature f2 averaged over a constant energy orbili^: 

-, , I dT dV 

(27r)^i^(£o) a|eo| dS 

Therefore, we obtain 

S{e,)ll^S{e)ll + \e,\^ (33) 
which does not require the cyclotron orbit to be circular. The energy quantization condition can now be rewritten as 

S{e)ll = 2^[n +\]- ^^Ml ^ 2n[n +\]- nWc (34) 
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in which we recognized the topological Berry phase. Inverting this last relation S{e)l'^ = 27r[n + (1 — Wc)/2] allows 
one to obtain the (semiclassical) Landau levels for the whole energy band. Finally, the Landau index shift is 

7L = 2 - — (35) 

and the winding number only matters modulo 2. This last equation is the central result of the paper. It shows 
that the Landau index shift 7^ is related to the topological part of the Berry phase nWc and not to the complete 
Berry phase r(C). The important point in the proof is the cancellation in the phase S{e)lg between the non- 
topological part of the Berry phase T{C) — ttWc ~ — |eo|(ir/(i|£o| and the orbital magnetic moment contribution 
AiBd{Sl^)/d\eo\ = \eo\dT / d\eo\. Physically, the topological Berry phase ttWc is just the usual tt phase that a bi- 
spinor acquires in Hilbert space as a result of a 27r rotation in position space. Here the spin 1/2 is actually the 
sublattice pseudo-spin. 

In the following, we consider several concrete examples such as boron nitride, graphene mono- and bilayer. These 
examples are treated either in discrete lattice models or in their continuum limit (effective low energy models). 



IV. EXAMPLE 1: TIGHT-BINDING MODEL OF BORON NITRIDE 



In this section, we consider a single layer of boron nitride, which has a honeycomb lattice with two crystallographi- 
cally and energetically inequivalent atoms (boron and nitride, usually called A and B) as a basis. Because of the two 
different on-site energies Ea — £b = 2A 7^ 0, the inversion symmetry is explicitly broken leading to a gap opening. We 
use a tight binding model, with hopping amplitude t and nearest-neighbour distance a, given by the following 2x2 
Hamiltonian in {A, B) subspace: 

H{k) = ^^^^ ^ with/(fc) = -t[e-''"-^i + 6-''=-^^ + e-''^-^-^] (36) 

where k is the wavevector in the entire Brillouin zone [A; = corresponds to the center of the BZ, i.e. F point], 
^i;^2,<53 are vectors connecting an A atom with its three nearest B neighbours and ai,a2 span the Bravais lattice 
[we follow the notations of Bena and Montambaux— : our H{k) corresponds to what they call basis 11^]. Note that, 
contrary to |/(A;)|, f{k) does not have the periodicity of the reciprocal lattice but satisfies f{k + G) = f{k) exp{iG-d3) 
where ^3 is the vector relating the two atoms A, B of the basis. This case exactly corresponds to that of section III 
with a specific form for f{k). The quantities of interest (Berry curvature, orbital magnetic moment. Berry phase, 
winding number) can be directly computed from the results obtained there. 




FIG. 1: Berry curvature Q [in units of a^] in the conduction band of boron nitride as a function of the Bloch wavevector 
{kx, ky) [in units of 1/a] in the entire Brillouin zone for A/t — 0.1. Tire lattice vectors have been taken as ai = ^ae^ + \aey, 
a2 = —^ae.x H- fae^. Left: three dimensional plot (fc^;, fc^, fi). Right: contours of iso-curvature in the Brillouin zone. 



The curvature is given by 



^ik) = j^;^ sin(fe . ^^) sin(fc . ^^) sin(fc • ^^) (37) 
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where |eo(fc)P = + |/(fc)p, see fig. 1. Note that the curvature has both the C3 symmetry and the translational 
symmetry + G) = il{k)) of the triangular Bravais lattice. 

The orbital magnetic moment is easily obtained from Ai ~ eeoil and is shown in fig. 2. 




FIG. 2: Orbital magnetic moment M [in units of eta^/h] in the conduction band of boron nitride as a function of the Bloch 
wavevector {kx, ky) [in units of 1/a] in the entire Brillouin zone for A/t = 0.1. Left: three dimensional plot {kx, ky,M). Right: 
contours of iso-A^ in the Brillouin zone. 

Because of time reversal symmetry, the curvature satisfies r2(— fc) = —^l{k) and its integral over the entire BZ 
vanishes. As inversion symmetry is absent fc) 7^ 51(fc). 

The Berry phase for a cyclotron orbit C of constant energy Eq is r(C) = TrlVcil — j^] where Wc = —a §^ (10/27: is 
the winding number, which is ±1 when encircling a valley (because of a vortex in 9) and when the orbit is around 
the r point, see fig. 3. A saddle point in the energy dispersion at jsol ~ vA^M^ separates the cyclotron orbits 




FIG. 3: Left: Isoenergy hues (eo(fc)=constant) of boron nitride in the first Brillouin zone for A = 0.1 [energies in units oit\. In 
the semiclassical limit, cyclotron orbits in reciprocal space follow the isoenergy lines. Right: Polar angle on the Bloch sphere 
6{k) = — Arg/(fe) in the BZ. The winding number Wc measures the topological charge of vortices in the polar angle 8. 

which encircle the two valleys from the cyclotron orbit which encircle the F point in the BZ. As a consequence, 

F(C) ^ -~a^n[l - A/|£o|] if A < |eo| < ^A^ + (i.e. Wc = -a£_ = ±1) 

= if v/A2 + <2 < < ^A2 + (3t)2 (i.e. Wc = 0) (38) 

We checked this simple expression for the Berry phase along a cyclotron orbit numerically by directly computing the 
integral of the curvature in k space over the area encircled by the cyclotron orbit. 



9 



Prom the energy quantization relation S{e)l'^ = 27r[n + 1/2] — ttWc it is now possible to obtain the (semiclassical) 
Landau levels for the whole energy band of boron nitride. It shows that the Landau index shift 7/, = l/2± 1/2 = 
(modulo 1) vanishes for cyclotron orbits encircling a single valley {K or K'). Whereas for orbits around the F point, 
it is 7L = 1/2 + = 1/2. 



V. EXAMPLE 2: LOW ENERGY MODEL OF BORON NITRIDE (MASSIVE 2D DIRAC FERMIONS) 

We now take the continuum limit of a single layer of boron nitride. The low energy effective theory close to two 
incquivalcnt corners of the Brillouin zone (called valleys K and K') is now given by a massive 2D Dirac Hamiltoniai>i^: 



cos /3 ^ sin /3e 



i/,(fc) = efc.., + A.. = |eo|( -cos/? ) ^39) 



where cos/3 = A/|eo|, sin^ = fc/|eo| with |eo| = VA^ + P > and < /3 < tt/2, e{k) = Arg{k^ + iky) and 
f^{k) = ^|fc|e^'^^'^'^-'. The wavevector k is now defined from the K or K' points and not in the entire BZ. The Pauli 
operator vector is defined as cr^ = {ax,£,(Jy) where ^ = ±1 is the valley index (■C = 1 corresponding to the K valley). 
The Permi velocity v = 3ta/2 has been taken to 1. The most general single valued eigenfunction with eigenenergy 
eo = a;|eo| is ip{r) = Uk{r)e^'^'^ where the Bloch spinor is 

, , / cos(/3/2) \ 

1"^-^ ^ Usin(/3/2)e»«« 1 ifa-+l 



-$sin(/3/2)e-'«'' 
cos(/3/2) 



if a = -1 (40) 



The Berry connection is given by: 

A = -a^sm^{/3/2)V0 (41) 
Upon integration over the circular cyclotron orbit C of radius k, we obtain the Berry phase; 

T{k) = -aC27rsin2(/J/2) = -a^7r(l - cos^) (42) 

and the topological Berry phase: 

ttWc = -a^TT (43) 



The connection can be rewritten as: 



Note that 



E = (45) 

which is a manifestation of time-reversal cancellation. The Berry phase depends on the magnitude of the gap, which 
means that it now depends on k and therefore on the magnetic field B. Two limits of interest are the "ultra-relativistic" 
limit (A/fc -^0, (3^ n/2) 

P(fc > A) w -<7r = ttWc (46) 
and the "non-relativistic" limit (/? sa fc/A — )■ 0): 

P(fc < A) w r(0) (47) 
when A 7^ 0. The corresponding Berry curvature is: 

1 dP A 
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It docs not contain a singular term, except when A 0+: $7 = when k and J7 —a£,oo when fc = 0. Details 
of the calculation are given in Appendix |D] 
The orbital magnetic moment is^ 



M = eson = (49) 

ZEq 



As a side remark, we note that this orbital magnetic moment leads to a valley-Zeeman effect in presence of a magnetic 
field. In particular, at the bottom of the band fc — ^ 0, the orbital magnetic moment is A^(0) = — ^e/2A and the valley- 
Zeeman gap would be 2Ayz = 2M{0)B. Some effects related to this valley magnetic moment are discussed in Ref. [22. 
Here, we would like to point out one more effect, which could be relevant for graphene in the quantum Hall regime. 
In graphene - which is gapless A = in the absence of a magnetic field ~ it is possible to imagine a self-consistent 
mechanism at finite B leading to a valley-dependent gap opening for the n = Landau level. Indeed asking that the 
gap leading to a valley magnetic moment is itself the valley-Zeeman gap A = Ayz leads to A = Hv/IbV^ (x \/B. 
This single electron mechanism is similar but not identical to that proposed by Lukyanchuk and Bratkovskjs^l, as can 
be seen from the different magnetic field dependence of the gap (square root versus linear). A valley splitting of the 
n = Landau level of graphene has indeed been observed in a strong magnetic field^. However it is not yet clear 
what is the relevant microscopic mechanism (for a review see Ref. [26l ). 

From the Onsager relation and the Berry phase just obtained, we find the energy quantization condition S{e)l% = 
27rn -I- 1/2 — nWc with Wc = ~a^- The area S{e) = n[s^ — A^] has the same functional form as S{C) = 7r[el — A^], 
but the two quantities differ by the term —2TTeoA4B. It is S{e) which is directly related to the Landau levels [and 
not S{C)]. By inverting ^(e), the semiclassical Landau levels are: 



= S-'[2neBin + 1 - 1^)] = a Ja^ + 2eB(n + ^ - ^) (50) 



The energy is therefore quantized as 



a 



VA2 + eB2n' (51) 



where n' = n -f (1 -|- q;^)/2 is an integer. This result agrees with the exact expression for the Landau levels ([7|), 
including n' = 0. Indeed, n' = n = implies a = — ^, which gives e = aA = — ^A. It is a bit surprising that a 
semiclassical calculation (including terms of order h) is able to recover exactly a fully quantum result. This is actually 
a peculiarity of massive Dirac fcrmions and docs not occur in more general cases. 



VI. EXAMPLE 3: LOW ENERGY MODEL OF GRAPHENE (MASSLESS 2D DIRAC FERMIONS) 

As another example, we consider the case of graphene, which is a two-dimensional honeycomb lattice of carbon 
atoms. It can be seen as the limit of boron nitride when the gap closes because the two carbon atoms in the unit cell 
have the same on-site energy. It is a zero-gap semiconductor and its low energy effective theory - close to K ov K' - 
is given by a massless 2D Dirac Hamiltonian: 

H^{k)^ivk-cT^^^vk[^^':^, g j (52) 

where 6 = Arg(fca: 4- iky) depends on the direction of the Bloch wavevector fe [here defined from the K or K' points] 
and ^ = ±1 is the valley index (C = 1 corresponding to the K valley). The 2x2 matrix is written in {A, B) space. In 
the following, we take the Fermi velocity v = 1. The most general single valued eigenfunction with eigenenergy ak is 
■0(r) = Uk{r)e^^''^ where the Bloch spinor is 

\uk,a) = 4s ( f J^e ) ffci = +l 



V2\i^ 
V2\ 1 



ff a = -1 (53) 



where a is the band index: a = +1 [resp. —1] corresponding to the conduction [rcsp. valence] band. The first [rcsp. 
second] component of the spinor is the amplitude on the A [resp. B] sublattice for both valleys and the area of the 
system was taken as unity. 
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The Berry connection is given by: 



A=-^Vue = -^ee (54) 



which shows that it is a pure gauge except for the singularity at the origin. Because of this vortex, it gives a topological 
(quantized) Berry phase: 

T = <j> dk- A = -a^n (55) 

which is independent of the cyclotron orbit. Here the winding number Wc = —a^ = ±1, where is the chirality of 
the massless electron. This allows one to rewrite the Berry connection as: 

A=^Vk0 (56) 

The corresponding Berry curvature is singular 

n^TS^ik) (57) 
and the Roth- Wilkinson relation ^ between the phase mismatch 7 and the Berry phase F gives: 

1 r l + < 

7=--- = ^.0mod. 1 (58) 



which is consistent with the Landau levels found by McClurc e„ = aV2neB. The Berry phase is non-zero here because 
of the band degeneracy (Dirac point) and despite the inversion symmetry being present (which results in = in 
the absence of band degeneracyi^) . A singular Berry phase in a system with inversion and time- reversal symmetry is 
a signature of the presence of a Dirac point. 
The orbital magnetic moment is also singular: 

M = er5{k)Si9) ^ -a(TTe5{k)6{e) (59) 

However, it plays no role in the quantization of cyclotron orbits for massless Dirac fermions because the area of the 
cyclotron orbit at constant energy iS'(e) = S{C) — 2TieoM.B is equal to S{C) as EqM oc k5{k) = 0. 



VII. EXAMPLE 4: LOW ENERGY MODEL OF A GAPPED GRAPHENE BILAYER 

The low energy effective theory close to K and K' of a gapped bilayer graphene is given by the following 
Hamiltonian^S: 

'll ) (60) 

where m is an effective mass and = Arg (k^ + iky). The function is therefore /j(fc) = (— A;^/2m) exp (i2^0), 



which shows that Oik) = 2^(/)(fc) - tt and |£o| = y + {k'^/2my. 

The Berry phase is r(C) = 7rM^c'[l — A/|eo|] where the winding number is: 

Wc = 2< (61) 

because the phase 9 rotates twice as fast as 0. From the previous analysis of the semiclassical quantization 
condition ([34]) including the effect of the orbital magnetic moment, we find the energy quantization condition 



S{e)l% 2TTn + 1/2 - nWc- Inverting S{e) = 2tt77iVP~~S^, we obtain the semiclassical Landau levels: 



a\lA^+ ujl{n + i - a£,Y ^ + u^{n' + i)2 (62) 

where the cyclotron pulsation lUc = eB/in and n' = n — is an integer. The quantum mechanical result^ is 
e„ = -|- uj'^n(n — 1), which agrees with the semiclassical result including the n*^ order. Here, however, the 

semiclassical results does not match the quantum result to all orders in h. 
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VIII. CONCLUSION 

We have studied wave effects in the semiclassical quantization of cyclotron orbits in coupled two-band models, 
focussing especially on the case of boron nitride. Two main results of the article are the following: 

First, although the phase mismatch 7(C) appearing in the Onsagcr quantization condition is related to the complete 
Berry phase r(C), the Landau index shift 71, only gets a contribution from the topological part of the Berry phase 
ttWc (winding number of the pseudo-spin 1/2). The latter is a topological invariant, which allows one to distinguish 
between two types of band insulators. On the one hand, zero topological Berry phase indicates that if inversion 
symmetry is restored, the bands are well separated and no Dirac points are present. On the other hand, a non-zero 
topological Berry phase is a signature of the presence of Dirac points in crystals with inversion symmetry. Therefore 
a shift in the Landau level index is related to a non-zero topological Berry phase, which signals the presence of 
underlying Dirac points (which are only revealed if inversion symmetry is restored) . 

Second, computing the Berry curvature in the entire Brillouin zone requires care in defining the fc-depcndent 
Hamiltonian. In particular this Hamiltonian should be written in what Bena and Montambaujs^^ call basis II and 
not in basis I, which is the basis that automatically emerges when performing the unitary transformation (|12p. The 
Berry curvature is a local physical quantity that in principle could be measured. A challenge would be to design a 
"Berrymeter" to measure this curvature in the entire Brillouin zone. An idea would be to measure the anomalous 
(/-factor, which is due to the orbital magnetic moment Ma and contains the same information as the Berry curvature. 
This could be done as a function of doping - e.g. electric doping in graphcne with a gate - giving access to local 
quantities. 
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Appendix A: Physical interpretation of the semiclassical quantization condition 

Here we give a physical interpretation of the semiclassical quantization condition for the cyclotron orbit. These 
results are certainly not new, but we collect them because they seem not to be so well-known. Physically Onsager's 
quantization is the condition for the single-valuedness of the semiclassical wavefunction. It states that the total 
stationary phase accumulated by an electron around its cyclotron orbit is the sum of four terms and should equal 
zero modulo 27r: 

(j}=: hkx 2TTr/n -eB X Trr^/ft + r(fc) - tt = 2Tm (Al) 

where n is an integer and — e < is the electron charge. These four terms are: the spatial de Broglic phase k x 27rr; 
the Aharonov-Bohm phase ~eB x irr'^/h; the Berry phase T{k); and the Maslov contribution of —tt. 

The two first terms are classical (they arrive at order hP). The de Broglie phase is just the accumulated phase 
of a (quasi) plane wave on a trajectory of length 27rr. For the cyclotron orbit, because classically hk = eBr as 
hk = —er X B, it can be rewritten as eB x 27rr^/^. The Aharonov-Bohm phase comes from the fact that the electron 
surrounds a region of non-zero magnetic flux $ = i? x tt?'^ and the minus sign comes from the negative electric charge 
of the electron. It is given by — 27r$/$o = —eB x irr^ /h, where $0 = h/e is the flux quantum. The Aharanov-Bohm 
phase can be seen as a Berry phase due to magnetic curvature in real space. Together these two terms form the 
classical reduced action (divided by h): Ad/h = § dr ■ p/fi = tt/c^ x h/eB = S{k)lg where p k, hk — eA. 

The two other terms are the first quantum corrections to the classical action (they appear at order h and represent 
wave effects): the Berry phase and the Maslov contribution. The Berry phase is due to curvature in fc-space because 
of the torus-like topology of the Brillouin zone. It can be seen as an Aharonov-Bohm phase due to a "magnetic 
field in k space" Jl, whose fiux is 17 x Trfc^. The Maslov contribution comes from two caustics (Maslov index of 
2) on the cyclotron orbit, each contributing a factor —tt/2. The caustics represent singularities in the semiclassical 
wavefunction, where the probability density diverges and the phase picks an extra —tt/2 factor. The caustics are 
actually not properties of a single orbit but of a family of classical orbits. For a detailed discussion of caustics and 
the extra tt/2 phase (in the context of optics) see Ref. [2^. 

Collecting these four terms, equation (|A1[) can be rewritten as: 

^=5(fc)/| = 2.[n+i-ffl] (A2) 
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which is precisely the semiclassical quantization of a cyclotron orbit including terms of order H, with 7 = 1 /2— r(fc)/27r, 
see Eq. ([TH)) . 



Appendix B: Hamiltonians 

In this appendix, we discuss the relation between the three kind of Hamiltonians used in the main text, namely H, 
H{k) and H{k). Hats are used to distinguish operators acting in the complete Hilbert space from those solely acting 
on band indices. 

1) The original Hamiltonian of the system is called H. In its eigcnbasis of Bloch states it reads: 

H = J2^"^o{k)\k,a){k,a\ (Bl) 

k,a 

2) Next, we define the unitary operator U{k) = exp(— ife • r) - where r is the complete position operator and not, 
for example, merely the Bravais lattice position operator -, which transforms Bloch states \k,a) into their it-part: 
U{k)\k , a) — \uk,a)- This transformation is just a translation by k in reciprocal space. Performing this unitary 
transform on H, we define the fc-dcpendcnt Hamiltonian: 

H{k) = U{k)HU{ky = ^e,,o(fc')exp(^(fc' - k) ■ f )|Mfc',„)(ufc',„| exp(-^(fc' - k) ■ r) (B2) 

k' ,a 

It is still an operator in the complete Hilbert space but it depends on A; as a parameter. This transformation actually 
defines a whole family of Hamiltonians (one for each wavevector k in the Brillouin zone). 

3) The 2x2 fc-dcpendcnt Hamiltonian is defined as the restriction of H{k) on the fixed k subspace: 

H{k) = P{k)H{k)P{k) = Y,^c.fi{k)\uk^o.){ukA (B3) 

a 

where P{k) = \ukM){uk,a\ are projectors on the k subspace. Note that H{k) is only an operator in band index 
subspace. In the case of only two bands, it is therefore a 2 x 2 matrix. Note that H{k) is not periodic in reciprocal 
lattice vectors but its eigenvalues are. 



Appendix C: Orbital magnetic moment for an electron- hole symmetric two-band Hamiltonian 

In this Appendix, we prove that there is a simple relation between the orbital magnetic moment and the Berry 
curvature in the case of a two-band model with particle-hole symmetrj*22. The 2x2 Hamiltonian is 

H{k) =Y,£oAk)Pc.{k) (CI) 

a 

where Pa{k) = \uk.a){uk.a\ arc projectors on each of the two bands (labeled by a = ±1) and eo,Q(fc) arc the band 
energies. Particle-hole symmetry together with time-reversal symmetry implies that 

£0-a{-k) = -£o,a(fc) = -£0,a{-k) (C2) 

and the Hamiltonian becomes 

iJ(fc) = £o,+ (fc)[P+(fc) - P_(fc)] = eo,-[P-(fe) - P+(fc)] (C3) 
Using the unit operator in the reduced fc-spacc I{k) = P+{k) + P^{k), we can write: 

eo,a(fc)/(fc) - H{k) = 2eoAk)[I{k) - Pa{k)] (C4) 
Therefore the orbital magnetic moment (in the upper band, e.g., a = -1-1) is 

M+ = t^{Vu+\x[eo,+ ~H{k)]\Vu+) 

= *^£o.+ (Vu+| X [/ - P+]\Vu+) (C5) 

But (Vu+I X /| Vu+) = — zrj+ by definition of the Berry curvature and (Vu+| x P+| Vu+) = (Vit+|u+) x (u+| Vm+) = 
,4+ X A.+ = by definition of the Berry connection, therefore 

M.a = T^O.Ma (C6) 
ft 
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Appendix D: Berry curvature of a massive Dirac fermion 



The Berry connection is given by: 



A^^-^VuO (Dl) 



Therefore: 



n = Vx A=^^fflx V0 + ^^VxV0 (D2) 
27r 27r ^ ' 



Using that V x V6' = 2TT5^{k)ez, we obtain: 

Therefore, if A 7^ 0, r(0) = and the Berry curvature is = (dr/dfc)/(27rfc) = -<A/(2|£oP). But if A = 0, 
r(fc) = r = -a^TT = r(0) ^ O, dr/dk ^ O and the curvature is 17 = r6'^{k). 
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More generally, we consider the case where a gap is opened in graphene's band structure as a result of inversion symmetry 
breaking. Two-dimensional boron nitride is the most obvious example. However it has a very large gap of roughly 2A ~ 6 eV, 
that is even larger than that of silicon dioxide, and it cannot be used as a semiconducto]>ii. A more interesting system 
for electronic transport is graphene on a commensurate (bulk) boron nitride substrate. It is predicted^^ to have a gap 
2A ~ 50 meV. In the following, we will assume that the gap is smaller than the bandwidth. 
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This result is valid for n > 0. When n = 0, the result requires the inclusion of both valleys at the same time: en=o ~ — C^, 
where ^ = ±1 is the valley index {K or K'). 

When studying an electron in a periodic potential in the presence of a magnetic field, one should be careful in defining 
various momenta. Here we consider four such momenta. First, there is the canonical (or linear) momentum p which is 
canonically conjugated to the position r. Second, there is the crystal momentum hq (where q is the Bloch wavevector) 
defined by the Bloch theorem. Third, in the presence of a vector potential, there is the gauge invariant momentum II which 
is obtained from the canonical momentum by minimal coupling II = p + eA (in the absence of a periodic potential, it is 
directly related to the velocity II = mv). Fourth, there is the gauge invariant crystal momentum hk, which is related to 
the crystal momentum by hk = hq + eA (in the absence of a magnetic field, the two are equal k — q and we therefore 
usually use fc). The gauge invariant crystal momentum is the one appearing in the semiclassical equations of motion. It is 
sometimes a valid approximation to neglect the difference between linear momentum and crystal momentum (think of the 
Peierls substitution), in such a case p ~ ftq = ftfe — eA. 

Note that in addition to the self-rotation A4q (fee) the magnetic moment of the electron also gets a more familiar contribution 
from its center of mass motion, which is taken care of by ea,o(fec ~ Qc + eA{rc)) ~ Safiiq^) + eA{rc) ■ Vfe^e^.o ~ 
^afiiq^) + ^(^c X Vfe^eQ,o) • B in the symmetric gauge where A{rc) = B x rc/2. The corresponding magnetic moment is 
— ^Vc X Vfe^£a,o ~ — f'"c X fc as expectcd. 

When the dispersion relation eo{k) is not isotropic, the cyclotron orbit in k space is not circular and the Berry curvature 
explicitly depends on the wavevector. Hence the necessity of defining an averaged Berry curvature. Another expression for 
this quantity is n(eo) ~ [(27r)^z^(eo)l~^c(r/d|£o| where i^(eo) = (27r)^^d5'/d|£o| is the density of states per unit area. 
When computing the Berry curvature and related quantities, there is no "choice of basis" for H{k) in the sense of Ref. [2^ . 
Indeed, the basis is fixed by the unitary transform exp(— ifc • f) used to define H{k) and this gives basis II. Writing H{k) 
in basis I, instead of II, amounts to replace /(fc) by //(fc) ~ —t[l + e^^^'""^ + e^**" "^]. To check that this is not correct, 
we computed the Berry curvature by boldly replacing / by // in the corresponding formulas and found a different Berry 
curvature, which did not have the C3 symmetry. 



